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' A difference equation w.r.t. space-time dimension d for n-point one-loop integrals 

. , with arbitrary momenta and masses is introduced and a solution presented. The result 

can in general be written as multiple hypergeometric series with ratios of different 
Gram determinants as expansion variables. Detailed considerations for 2— ,3— and 
4— point functions are given. For the 2— point function we reproduce a known result 
in terms of the Gauss hypergeometric function 2^1- For the 3— point function an 
expression in terms of 2^1 and the Appell hypergeometric function Fi is given. For 
the 4— point function a new representation in terms of 2F1, Fi and the Lauricella- 
Saran functions Fs is obtained. For arbitrary d = 4 — 2e, momenta and masses the 
2— ,3— and 4— point functions admit a simple one- fold integral representation. This 
representation will be useful for the calculation of contributions from the e— expansion 
needed in higher orders of perturbation theory. Physically interesting examples of 3— 
and 4— point functions occurring in Bhabha scattering are investigated. 
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1 Introduction 



The calculation of radiative corrections in the electroweak Standard Model (SM) is espe- 
cially demanding because of the many species of particles and fields and the large variety of 
interactions between them as well as the many mass and energy scales which typically show 
up in a high energy scattering process. Practical problems in the calculation and numerical 
evaluation of Feynman integrals were encountered at an early stage of SM calculations be- 
yond the tree level, at the level of the one-loop integrals already. To a large extent, known 
analytic results and techniques have been reviewed, extended and discussed long time ago 
in PP . In addition M. Veltman first had the idea to develop a library of numerical routines, 
the program FormF |2|, which allowed to calculate the one-loop radiative corrections for 
2 — >■ 2 scattering processes in a numerically stable way. Unfortunately, this program was 
only running in CDC assembler, and used very high precision (more than 100 digits) to 
overcome the numerical instabilities that plagues the calculation of one-loop electroweak 
radiative correction. Later G.J. van Oldenborgh implemented the most relevant one-loop 
integrals as a normal Fortran library FF |T, which utilizes alternative algorithms which 
run in double precision Fortran 77. This program was very useful for many of the calcu- 
lations which were needed for precision physics at LEP. Meanwhile a number of extensions 
to two- and more-loops have been developed. For example, the programs SHELL2 jl] and 
MINCER allow to calculate special types of integrals. 

More recently we face new problems in view of the requirements for future colliders. For 
example, we certainly will need full electroweak two-loop calculations for 2 — > 2 fermion 
processes as well as full one-loop 2 — > 4 fermion processes for precision physics at a future 
e^e~ linear collider like TESLA. In both cases the existing program libraries need important 
extensions. 

In the first case, for the renormalization of two and more loops, one needs one-loop 
integrals to higher order in the 5-expansion (rf = 4 — 2e), while in the second case efficient 
algorithms are needed for the calculation of 5—, 6— and higher- point functions, which at 
least on the level of the scalar integrals are ultraviolet finite in c? = 4 dimensions. Usually, 
one may reduce them to a sum of 4— and lower-point functions, however with "unphysical" 
external kinematics. While FF serves all requirements to calculate 2 — 2 processes at one- 
loop, it does not work in general if we try to evaluate the expressions we obtain as a 
result from the reduction of higher point functions and we thus have to extend this tool 
appropriately. 

Another important extension we need stems from the fact that unstable particles, like 
the massive gauge bosons W , Z, the Higgs and the top quark, show up as intermediate 
states in scattering matrix elements. Most interesting is their production near resonance, 
where ordinary perturbation theory breaks down and requires to work with finite width 
propagators not only for single resonant virtual particle lines but also within loops. Most 
obviously this shows up in the interplay between virtual and real photon emission in infrared 
sensitive quantities. 

In this paper we propose a new method to deal with one-loop integrals which allows for 
extensions in several respects. The algorithm we propose works for arbitrary space-time 
dimensions, for complex masses and momenta and allows to find the appropriate physical 
domain of the analytic structure. 
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The paper is organized as follows: in Sect. 2 we recall the basic recursion relation and 
introduce the relevant kinematical quantities. Based on this recursion relation a difference 
equation and its solution is given in Sect. 3, which quite generally represents n— point func- 
tions in terms of series over (n — 1)— point functions. An asymptotic method to determine 
the 'boundary term' is discussed. In the following sections, we give explicit representations: 
in Sect. 4 for the 2— point function, where also the approach itself is demonstrated in great 
detail; similarly in Sect. 5 for the 3— point function and in Sect. 6 for the 4— point function. 
Explicit examples needed in Bhabha scattering are also presented. For the boundary term 
of a 4— point Bhabha diagram a differential equation is set up as alternative method for its 
determination. 



2 General remarks on one— loop functions 

Relations between Feynman integrals in different dimensions are known since quite a long 
time [0] . Recurrence relations for n-point one-loop integrals, required for our investigation, 
are given in [7|, simplifying and extending the work of .8J. We consider scalar one-loop 
integrals depending on n — 1 external momenta: 



d'^q " 1 
11 



3=1 S- 



(1) 



where 

Cj = {q - pjf - m] + ie. (2) 
The corresponding diagram and the convention for the momenta are given in Fig.l. 




Fig. 1: One-loop diagram with n external legs. 

A recurrence relation connecting integrals with different space-time dimensions was 
given in [7j 



i=i 



2A„ + 5;](afcAOk- 



k=l 



Ad) 



(3) 



where dj = d/dm^, and 
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A„ = A„({pi, mi}, . . . {pn, rrin}) 



Yii Yi2 . . . Yin 
Yi2 Y22 ■ ■ ■ 

Yin Y2n ■ ■ ■ Yjin 



(4) 



(5) 



Gn-l = Gn-l{pi, ...,Pn) 



{Pl-Pn)iPl -Pn) 
{Pl - Pn){P2 - Pn) 



{Pl ~ Pn){P2 - Pn) 
{P2 - Pn){P2 - Pn) 



{Pl -Pn){Pn-l -Pn) 
{P2-Pn){Pn-l -Pn) 



, (6) 



(Pl - Pn) (Pn-1 - Pn) (P2 " Pn) (Pn-l - Pn) ■■■ {Pn-1 " Pn) {Pn-1 " Pn) 

Go = —2, 

Pi arc combinations of external momenta flowing through i-th Hnes, respectively, and 
is the mass of the i-th line. In the following we also quite often use the abbreviation 
Plj — (Pi ~Pj)^- By shifting the integration momentum we choose in general Pn — 0. Where 
no confusion can arise, we simply refer to the above functions as A„, Gn-i- Considering 
integrals for n = 2, 3, 4 in the next sections we will use also an indexed notation for A„ and 

Gn-l 

Ki2...in = ^n({Ph,mil}, {Pi2,m2}, {Pin^^n})^ 



9iii2---in Gn—liPi\iPi2i ■ ■ ■ iPin)' 



(7) 



The indexed notation will be useful when considering integrals obtained from /^^^^ by con- 
tracting some lines. Rather frequently the results will depend on ratios of Aj^jj-in 



9; 



and therefore it is convenient to introduce the notation 



With this definition a useful relation is: 



^ij...k 



9ij---k 



•^fc — l^fc + 1 •••^71 

Using 

n 

i=i 

one shows that to all orders in e 



(8) 



(9) 



(10) 



^ni2...jn(W - ^e}) = Ki2...in{{'l^l}) + Whi2...in 



(11) 
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and therefore the e prescription for r is rather simple (with the same e for all masses) 



J J 



ie. (12) 



As we will see later, the results for one loop integrals will be expressed in terms of hy- 
pergeometric functions depending on ratios of different r's plus terms proportional to j. 
Relation (fT^ makes the r's appear like masses, i.e. their ie prescription is the same as for 
masses (their dimension is mass squared). This simple property is of extreme importance 
for the study of the analytic behavior of one loop integrals, in particular their analytic con- 
tinuation to different kinematical regions for arbitrary space-time dimension. For brevity, 
we shall omit "causal" ie's below. 

3 Solution of the difference equation 

In this section we will derive an explicit solution of relation Q for z/j = 1. This relation 
has the simpler form 



{d-n + l)G„_i/(^+2) 



2A„ + ^(afcA„)k- 



k=l 



(13) 



If we assume that, evaluating n -point functions, we already know n — 1 point functions, then 
relation ()13|1 represents an inhomogeneous first order difference equation with respect to d. 
Methods to solve this kind of equations are well described in the mathematical literature 
0, [ID]- By the redefinition 

" r(^) [Cn-J " ^ ' 

we obtain the simpler equation 

p / d-n+l \ „ \ - n 

= + -L-L-^ %1 ' E(*A..)k-^'. (15) 

Without loss of generality we can parameterize d as 

d = 2l~ 2e, (16) 

where / is integer and e small and possibly complex. Then the solution of the equation for 
I^^^ can be written as 

I r (r — ] ~ F — /n \ r-l-£ n 

= — — ^ %^ i:(aA.jk-e--->+i;„(^), m 

r=Q ■^'-^n \ '-^n ' k=l 

where 6„(e) = In'^ ^"^"^ is an / independent constant. With J2l=o = J2'^o~^r^i+i 
shifting the summation index in the second sum like r — » r + / + 1, the solution ()17p can 
be rewritten in the form 

oo r (r I d-n+l \ „ X r+l n 

If-"' = - E s A i:(aA„)k-e«) + 5„(.) (18) 

r=0 ^ ^ k=l 



by redefining the / independent 'boundary term' The final result for I^^ then reads 

As we will show in the next sections, 6„ can be determined from the asymptotic behavior 
of for d ^ cxo or by setting up a differential equation for it. This term depends on the 
kinematic domain. 

The series in the above solution converge in general if the expansion parameter does not 
exceed 1. If it does, one has to continue the series analytically. This can be done by different 
methods. If the result is already obtained in terms of known hypergeometric functions, e.g., 
then known formulae for their analytic continuation can be applied. Another method will 
be to modify the procedure of solving the difference equation. According to the general 
theory of difference equations, we should repeat the derivation in our case by using the 
following parameterization for d (changing the sign of /): 

= -2/ - 2e. (20) 
With this parameterization for the function 



= (21) 



we obtain from Eq. 



Y ( -2l~2e-n+l \ „ . -l-e-1 n 

= fi- ^ ^) Ud.A^)^-lL-''-'-'''- 



V A„ / 



(22) 



Solving this equation we obtain for I^^ 

- M.) - jj^^ t (^) LA^r\-it-, ,23) 



where {a)r = T{r + a)/r{a) is the Pochhammer symbol. Here we have for convenience 
introduced the same notation for the boundary term as before, keeping in mind, however, 
that the value may be different from the one in ()19|) . In order to obtain convergent series 
for all the contributions in J^k^ can use different parameterizations for d in the various 
terms. 

In general the situation is such that to obtain the multiple series is straightforward. The 
problem is rather to determine the boundary term. In the next chapter we discuss for this 
purpose the asymptotic behaviour for d oo. 

3.1 Asymptotic behavior of I^^ for d ^ cx) 

For arbitrary n the asymptotic value of I^^ as ^ oo can be found by utilizing the 
parametric representation for i.e., with 

dXi . . . dXn-l r(''^) X2 ^ ■ ■ ■ Xn-2 (24) 



ciC2...c„ Jo '"Jo [cixi . . C2X1 . . . x„_2(l-x„_i) + . . .+c„(l-xi)]"' 
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shifting q in order to remove the hnear term and integrating over q by means of 



d'^a 1 T (a- 

' 1) ^. I (25) 



[i7:^n]{q^-mir ' ' r(a)(m: 
any n-point function can be represented as a multiple parametric integral of the form: 

l'f^=V (n-^l / .../^ (ixi...rfx„_i/(xi,...,a:„_2)(/i„({pj,m,},xi,...,x„_i))^"". (26) 
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In the analyticity domain of the polynomial ms}, Xi, . . ., x„_i) > and 

is an integral of Laplace type. In the domain of nonanalyticity there are subdomains of 
the integration region where hn{{pj, ms}, xi, . . ., x„_i) < and due to that the integral I^^ 
gets an imaginary part. For this kinematical configuration I^^ may be represented as 

/W= r (n- y{dx}e{K)f{{x}){K{{p,,ms},{x}))^- 

+ cos^(rf-2n)^ / _{dx}f{Xi,...,Xn-2) \K{{pj,ms},Xi, . . . ,Xn-l)\^~'' 

±isin ^(d - 2n) ^ / {dx} f {xi, . . . , Xn-2) |/i„({pj, mj, xi, . . . , x^-i) 1 > .(27) 
^ 3 •'^j J 

Here VL~ are subdomains of the integration region where hn < 0. The sign of the imaginary 
part has to be determined from the ie prescription in hn- 

Our parametric integrals are of multiple Laplace type which in general can be written 

as 

F{\)= [ f{x)exp[\S{x)]dx (28) 

(identifying S{x) = ln(\h{x)\) and A = | — n ), where f2 is a bounded simply connected 
domain in k dimensional Euclidean space, x = {xi,...Xk) and S{x) is a real function. 
The functions S{x) and f{x) are sufficiently differentiable functions of their arguments 
throughout Q. Supposing that a relative maximum of S{x) in Q is achieved at the interior 
point X = x , ^ then at A — » cxd [TT] 

F(A) ~ exp[A^(lr)] ^ arX'"'^ (29) 

r=0 

(correspondingly this holds for a relative minimum). This expansion may be differentiated 
w.r.t. A any number of times. The leading term of the expansion is 

F(A) = exp[A5(x)] (27r/A)t /(|g) + ^(^^ ') ^ (30) 

|detS'^.^(x)| 



''if the function h„ reaches its maximum value on the boundary of the integration region then 
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where detS'^j^: is the Hessian determinant, the indices appearing as second derivatives. The 
behavior of I^^ for — > oo can thus be obtained. An extremum is found from 

dh„ , . 

and using the determinant representation of h^, it is possible to show that ^ 



and at this point 

hni{pj,ms},Xi, . . .,Xn-l) =ri,„ri- (33) 
All Xi are inside the interior region, i.e. 

0<^i<l, (i = l,...,n- 1) (34) 



if all derivatives 9fcA„ have the same sign. Conversely, it is easy to see that, provided is 
true, due to (fTUj) all 5fcA„ have the opposite sign as Gn-i- assuming Ylk=i ^k^n > 0(< 0), 
from the left hand side of the inequalities (jMj) follows J2k=i dk^n > 0(< 0), i = 1, . . . , n — 1. 
Multiplying ()H4|1 with the positive (negative) denominators of x,, the right hand side yields 

dn-i+l^n > 0(< 0),Z = 

The general idea to determine the boundary term 6„ from the asymptotic behavior as 
0? — > oo is as follows: from (fT^ and (fTHj) we see that for large d 

bn ~ rf^.^n- (35) 

Such a contribution can come from the asymptotic behavior of I^^ on the l.h.s. of (|T!I|l due 
to ()29|1 and This may happen, however, only if all < irj < 1, i.e. if for the specific 

kinematics under consideration an extremum occurs inside the integration region and an 
absolute maximum of does not occur on the border. Finally we can write in the region 
where the integrand of ijf^ has an extremum inside the integration region 



JW = -(27r) 



ri 




iri 




d 

1 rn 


' — n — 1 
2 




■p j d—n 
V 2 


^) 








(36) 



where r„ = ri ^. To determine the boundary term, convergent series are requested on the 
r.h.s. of (jHUj) . These series are convergent in general only in certain kinematical domains 
and only for these is the obtained boundary term the correct one. The method of analytic 
continuation to other domains has been indicated above. In fact, the summation with 



^There are other solutions, but for them in general one of the Xi{i = 1, . . .n — 2) equals such that in 
these cases no contribution to J^'^'' is obtained 

^If the function /i„ reaches its absolute maximum value on the boundary of the integration region then 
I^"^ has no asymptotic behavior like ~ and = 0. 
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I > and / < 0, respectively, for the analytic continuation turns out to be easier in general 
than the direct analytic continuation of hypergeometric functions. In the next section we 
demonstrate the procedure in all details for the 2— point function and also give an example 
for the summation with / < 0. 

Analytic continuation of the (generalized) hypergeometric functions can also be useful 
to find the asymptotic behavior, like double logarithms, for certain diagrams in specific 
kinematic regions. 



4 2-point function 

We find it convenient to label the lines of 12'^^ as i,j because the integrals 12"^^ will be 
encountered in calculating I^^^ (and I^^^) as a result of contraction of different lines. Ex- 
pression (jl9|) for l2^^ includes two one-fold sums over tadpole integrals given in (j25p . 
Ii'''' can also be obtained from (|19p assuming that in dimensional regularization Iq^^ = 



-r 1 



(37) 



For n = 2, substituting ()37|1 into p9|l gives (observe the different normalization of 62 in 
comparison with (HH)): 



(d) 



+ 



+ 



diXij ( 2 



(m 



r=0 



(i). 

2 



Ail 



(38) 



where 

and the infinite series in 



can be represented as hypergeometric functions, i.e. 



(39) 



2^1 



d-1 . ' 

d . ^ 

2 ' 



(40) 



The parametric formula for l2^^ is a one-fold integral 



4'^) = r(^2-^^ J^'dx, \p^^xl~x,{p: 



2 2 I 2\ , 21 #-2 



(41) 



The extremum of ho is located at 



Xi 



di Xij 
diXii + 82 Xi 



pf - — + 



(42) 
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where 

= -{pI? -mt- + 2plml + 2plm] + 2m^^m]. (43) 

and 

d,\ij = 2{pI - ml + m]) , djXij = 2{pl + -m]). (44) 

Since 

a maximum of /12 inside the integration region exists for Euchdean momenta pj^ < —{rrij — 
mf) (without loss of generahty we assume here rrij > rrii). A minimum exists 'inside' for 
pfj > — 1^1) ■ In the first case we have diXij < and d2\j < while in the second 

case we have diXij > and d2Xij > 0, i.e. in both cases the two derivatives have the same 
sign. 



350- 
300- 
250 
200- 
150- 
100 
50 H 




Fig. 2: Function hf^ ^ at p| = -3, mf = 2,m'j = 1, d = 18. 

For completeness we also give the imaginary part of 12'^^ on the cut, which we obtain 
from (El: 



(46) 



Im =±r(2- fj sin ^^^^(-4.)(^^) dx,{{x, - x^){xt - x,))^- 



where 



are solutions of the equation 



X 



(47) 



h2{{pij,mi,mj},xi) = 0. 



(48) 
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For p^. > {rrii + rrij)'^ both < xf < 1. Performing the integration in (jlHI), we obtain 



Im L 



id) 



Ptrr {d - 2) 



d-3 



(49) 



in agreement with the resuh jl2j obtained by another method. Determining now 62, we 
have to investigate the asymptotic behavior of l2^^ for large d. For convenience we give the 
following expressions for r^: 



{mi + mjy 



Pi 



[mi 



nii 



(50) 



and 



Tij - rrij 



We begin with pf- < —{m] — 



{diXi 



{djXi 



(51) 



16p^ ' ^ 16pj= 

< 0. Here we have a maximum like the one shown in 



Fig. 2 with < 3;i < 1. In this case rij > mf^rn?- so that the hypergeometric series in 



d/2 



converge and a term 
we can write the result in the form 



~ Tjj comes from the asymptotic behavior of I2''' ■ Evaluating this, 



+ 



a A 



H \ " 



(i). 



d-2 
, 2 



diXij 



+ 



J 

I'ij 



d-1 
2 



(52) 



For 



(m^ 



< < we have again Vij > m^,m^, but xi < , so that there is no 



~ r, 



d/2 



comes 



extremum inside the integration region and as a consequence no contribution 
from I2 and the boundary term is zero. In fact in this case the two contributions in the 
boundary term of cancel and the formula is still valid. 

Investigating as next p?- > we can make the general observation that Vij < rrii, m3, but 



rij < is possible . Thus, continuing further to < • < (m 



-rrii)'^ we have —00 < Vij < 0, 
but independent of its absolute value this does not cause problems since even if the series in 
()52|) do not converge, for negative arguments the hypergeometric functions remain in their 
domain of analyticity (also the boundary term vanishes) and ()52|) is valid if written in the 
form 



r(i-i) "r(^)^-- 



+ 



djXij 



+ 



diXij 



1 2 

2 



2 ' 



m- 



+ 



2F1 



d-l 

d^. 

2 ' 



(53) 
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By means of ^3] 



1 ^Zli ■ ' 
' 2 ' 7 

2 ' 



1 t^-i 
i, 2 

3 . 

2 ' 



a compact expression is obtained from (|^: 



+ 



r (i) r (i 



2-d 
2; 2 



|arg(l — z)| < TT 



(54) 



\ij L 



id) 



diXij 



r(2-f 



(m2)i- 



"T 2 



d-1 

2 
3 . 

2 ' 



*j J 



d 2 
2 



1 ^ni ■ 

' 2 ' 1 
3 . 

2 ' 



(55) 



i.e. the series for the hypergeometric functions 
do still converge, but beyond that it is a valid representation also when we proceed to 



which applies as before for rjj > m^, 



rrii 



< pfj < {rrij + since < rij < mf, m'j in this case, i.e. the argument of the 



hypergeometric function gets negative and therefore is in the domain of analyticity. 



Above the threshold, i.e. 
— s> +00, behaves as 



{rrij + rriiY < pfj, we have again rij < 0. Due to for 
—00. Since in this case < xi < 1 and h2 has a minimum. 



apparently for some value of pjj we have \rij\ > m'^,m'j and the modulus of the argument 
of the hypergeometric function in (j53p is < 1. Therefore, according to our approach, the 
boundary term given in (j53|) is also obtained in this kinematical domain. 
A very useful transformation-formula, applied to ()53|) . is 



2F1 



d-l 

2 
d . 

2 ' 



iF, 



2 

d . 

2 ' 



(56) 



i.e. for all Vij < the arguments of the hypergeometric functions can be transformed into 
values of modulus less than 1 such that the series converge 

To conclude, we have found for all 
except for < Vij < mf, even convergent series. If we wish to obtain a convergent series 



-00 < pjj < +00 a valid representation for /- 



2*^^ and 



in the latter case as well, we have to proceed as described in Sect. 3, (j23j) . i.e. to perform 
the summation with / < 0. Without repeating the calculation in detail, we here just give 
the result: 



(d) 



r(2-f) 



diXij 



(d — 3)(m| 



r 



3-d 
2 



r 2 



2 



Fi 



diXi 



5-d . 2 

— , ml 



djXij 



mr — r,- 



djXij 



d-3 
, 2 



{d — 3){mf 



2F1 



1,2- f; r,,- 

5-d . 2 



(57) 



Using 



2F1 



1 2 — - • 

b-d . Z 

2 ' 



{d - 3) 2F1 



1 2 — - • 

3 . 

2 ' 



+ 



r {^) T (1 



d-3 
Z 2 



largfl — z)\ < n. 



(58) 
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yields from (j57j) the already known result |14j : 



9ij I; 



(d) 

2 



r(2-f 



djXij 

43T 2 



1 4— d . », 

3 . -L 2 

2 ' 



43T 2-ri 



fm?) 2 



1 4— d . », 

3 . ^ 2 

2' "^i. 



(59) 



Having thus obtained different representations for l2^^ by different ways of solving the 
difference equation (see Sect. 3), we can also show that they agree: by means of 



2F1 



1 ^Zli ■ 
' 2 '1—7 

2 ' 



1 I 

' 2 ' 1 — - 

3 . ^ 

2 ' ^. 



(60) 



dSni) can be obtained from 

For future applications we also give the e expansion of 2-^1 occurring in (see [T3] 
formula A. 3 in the Appendix): 



2F1 



1 i- 

d . 

2 ' 



(1 - 2£)i/ 



+ ln(l - y)e - (Lhiy) + \n\l - y)) e' 
1-y V / 



+ ( 2S,2iy) + 2 ln(l - y)U2iy) - Li^iy) + - -y))e' 



(61) 



where 



y 



X 



1 + Vl^ 

The e expansion of 2-^1 in ()58|) can be obtained from 



(62) 



12 — -- 

' 3. ^' 1-^ 
2 ' 



1 



z{2 - d) 



1 1 



2 ' 



More information on the e expansion of hypergeometric functions is given in ^H] • 



5 3-point function 



In complete analogy to ()38|) we proceed for the 3-point function. In this case we have to 
sum according to ()19|) 2-point functions, which in Sect. 4 were represented in terms of 
hypergeometric functions 2-^1- In order to simplify the summation, we eliminate d in the 
first argument of 2F1 by means of 



2F1 



a, b : 



[l-z)-%F, 



a,c — b; 

c; z 



(64) 



Dissolving the result into double series, it can be written as Appell hypergeometric functions 
Fsiid - 2)/2, 1, 1, 1/2, d/2; x, y) defined as 



m,n=Q 



(65) 
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which in turn can be reduced to Fi by means of ^3] 

Fsia, a', /3, (3\ a + a'; x, y) = {1 - yY^' i^a, (3, /?', a + a'; x, j , (66) 
Fi being defined as 

F,{a, f3, f3', 7; x,y) = £ («)-+»(/^)-(/^0n ^^^. (g^^ 

In the next step the boundary term in ()19p has to be determined by means of the asymptotic 
method described in Sect. 3.1 . The polynomial (see (|27|) for n = 3) needed to determine 
a possible extremum 'inside' the integration region, is given by 

= -0:1X2(1 - - ^1X2(1 - X2)p-j - Xi(l - Xi)(l - X2)p% 

+ XiX2m'^ + xi(l - X2)m^ + (1 - ^i)^k- (68) 



It has an extremum at 



{diXijk + djXijk) _ _ diXijk 



^1 = V ' ^2 = ^ , , (69) 



provided the determinant 



d^h d^hs ( d^h3 Y_ jdjXijk + djXijky ^ ^ .^q. 
9xi(9xi (9x2i9x2 1(9x1(9x2/ 2G2 



(71) 



Further, the following second derivatives determine the type of the extremum: 

dxidxi {diXijk + djXijkY ' 
d^h^ ^ 2pl{d,Xijk + djXjjkf 

9X2(9X2 Gl ' ^ ^ 

If both are > 0, has a local minimum, if both are < 0, has a local maximum. Thus 
we see that a maximum will be achieved if 

^2 < 0, rijk> ml,m],ml, < 0. (73) 

Here the explicit appearance of nik in (ffT|l and of pfj in (ff2|) are due to the particular choice 
of (jUHjl . i.e. the choice of the Feynman parameters. In general no preference of any choice 



of parameters can occur, hence the above condition for Vijk and p'^j. 
The result for I^^^ finally reads ^7] 

— 7 — "' ,i\ ^3^^ = ^3 + 0ijk dkXijk + dkij djXijk + djki diXijk-, (74) 

r(2-f) 

where 



3 I d-2 

63 = 25 7r ^-g^Jk Tijl , (75) 
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(see (j36j) ) provided an extremum of {G2 < 0) occurs inside the integration region of the 
Feynman parameters. Otherwise 63 = 0. For Xij 7^ we have 



^ij ^ijk 



OiXij ^ djXij 



L ■ 

1 i 



1 - 



4F 



2F1 



1 d—2 . /V, 

d-1 . 

2 ' ^ijk 



^ (m^) 2 djXjj / d-2 d mf mf 



2{d-2) rz 



2 
z_ 



2 2 r^jT^ r^j 



+ 



(76) 



To our knowledge there exists no simpler hypergeometric representation of the 3— point 
function for d dimensions in the literature. For another approach see e.g. ^Hl- The 
function Fi admits a simple integral representation: 



Fi(a,/3,/5',7;a;,2/) 



r(7) 



1 u'^-^il-uV' 



-a-l 



T{a)T{-i -a) Jo {I ~ xuf {I - yuf 



du. 



which in our case is 



'd-2 Id ' 



d-2 



d-4 
U 2 



Jo (1 — xu)^yl — yu 



du. 



(77) 



(78) 



For further information on Fi, we mention for its analytic continuation jTHI and concern- 
ing methods for its numerical evaluation ^20j. For Xij arbitrary, using the transformation 
formula 



Fi 



'd-2 Id ' 
,1, 2' 2'^'^, 



1 1^—2 . rp 

i, 2 1 

— ' y 



+ 



(2 - d) ^/r^ 



(1 - x)y 

the following compact form for Oijk is obtained 



Fi 1,2--,1,-;1-- 79 
V 2 2?/ y{l-x)J 



9ij ^ 



ijk 



imy) 2 



2 1 



rnk 



( d 3 



+ diXij- 



{mj) 2 



— Fi (l,2-ll,|;l 



2(1-^ 



where again Fi has a simple integral representation 



Tij fij TJlj 



(80) 



Fi h, 2-^,1, |a;,y 



1 (iu (1 — xu) 2 
2io v^r^ (1 - yu) 



(81) 
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To find /i'^'* up to first order in e 



(4 — (i)/2, three terms in tlie expansion of tliis function 




Li2(-A) 



\n{A) 



+-ln(A)ln2 



wfiere 



A 



1- i 

X 



B 



(82) 



(83) 



Tliis expansion will be of interest for an evaluation of one-loop counterterms inserted into 
vertex functions, required in two-loop calculations. 

The general result ()74|) for the 3-point function with arbitrary masses and external mo- 
menta is finally a bit lengthy. To demonstrate the usefulness of our result, we choose as 
a particular case the on-shell 3-point function occurring in Bhabha scattering: Jq'' with 



ml 



0, ml 



m and p| 



,2 

Am? 



m". For < pfj < 4m we have G2 < 0, i.e. an 



-7-^ — 2- > 1, such that it is not 'inside'. Therefore 



extremum (minimum) exists, but Xi 
63 = in this case. 

Evaluating the three remaining terms in (ff^ . we put = 1,2,3, i.e. we have to 

calculate 6^123,^312 and 6*231. With A12 = — P12 and A31 = A23 = we see that for the two 
latter cases we have to use (jHOl)- In the first case only the hypergeometric function 2F1 
remains: 



'123 



d-l 
2 



P12J 



2^1 



d-2 



' 2 
d_^l 

2 



1 - 



P?2 



-TT- ; 4m2 



(84) 



For 6*312 and 6*231, using (jHUj) . only one of the Fi functions remains, respectively, since 
c^sAsi = and 93A23 = 0. Moreover the two last contributions in ([7^1 are the same. 
Evaluating thus 6*312 we obtain: 



e 



^ 2^ ^ 4m2 



d 3 

^'^'^~2'2'^"4m2' 



P?2 



^2(^2 \2 Q^\^^^ 



With A 123 = -2m\pi^j 

Fi{a, b,b ,c, w, 1) 



-2pl2, <9iAi23 

rfc)r(c- 



C^2Al23 
-b' 



Am p 



^-"^ and 



125 



T(c-a)T(c~b'] 



2F1 



a, b ; 



w 



(85) 



(86) 
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the final result is 



J, 



G 



(m 2 



1,1; T 

d-l . -L 

2 ' 



r(2-f) 2{d-3) 
The result for Jg sX d = A was given in 



p2 




m 


(- 


2N ^ 


Am? 




2r 


m 


m? 



d-2 . 

2 ' 1 
d-l . -L 

2 ' 



4m2 
(87) 



6 4-point function 

The calculation of the 4— point function follows the same scheme as has been described in 
detail in the foregoing sections. Concerning the boundary term, it is in principle not difficult 
to find a relative extremum in a cube, but to find a proper formulation of all conditions (as 
has been done in the case of the 3— point function, see (fTnjl - (fTSjl ) in terms of 'standard' 
expressions like G3 etc. turns out to be tedious. However, due to the discussion after 
one excludes an extremum 'inside' if not all derivatives dk^n have the same sign. In the case 
of the 4— point Bhabha diagrams, nevertheless, we demonstrate an alternative approach to 
determine the boundary term. 

We label the lines of the box integral in natural manner as i,j, k,l . As one particular 
case we select an I^'^^ obtained from I^^'' by shrinking line / and take into account at first 
only the first Fi term in (fTBj) : 

/ diXijki\ ^ I d — 3 



Ad) 
'4,Fi 



E 



9ijkl 



j{d+2r) 
'3,Fi • 



(88) 



\ijkl ) .r=o \ ^ / r \ ^i'^' / 

In analogy to (jMj) for the 3— point function, to simplify the summation, we get rid of d in 
the first argument of F\ by means of the transformation formula pH] 



Fi(a,/5,/5',7;a;,y) = (l-x)-'^(l 
After this simplification we arrive at 



X 



y 



X — 1 y — I 



(89) 



8 \ A, 



ijkl 




^ijk 



E 



m+n 



^ijk 
(l).(l)m(i) 



rr 1- 



r,m,n=0 I 2 / 

1 ^Aijkl 




2 



^ijkl , 



2 

mf 



nit 



fijk , 



nit 



(mf 



d-2 
2 Tijk 



rr 1- 



mt 



^ ^ ^ \ d d d mf 
1 !■ 1 1 — ■ — — — ■ — — — 

2 2 2 2 Tijki mi - Tijk mi - Vi 



(90) 

zjkl '"'i I ijk '"'i '13/ 

where the Lauricella-Saran function Fs is given in terms of a triple hypergeometric series 

I22I, E 



Fs{ai, 02, "2, A, /52, 71, 71, 71; ^, y, z) 

{ai)r{a2)m+n{l3l)r{(32)m{l3'i)n 



E 

r,m,n=0 



(7i)r+m+n r\m\n\ 



x'y^'z''. 



(91) 
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In our case the following integral representation of Fs is useful 

r(ai)r(7i - tti) 



r(7i) 

1 ^7i-cn-l^-|^ _ 



/o {l-x + tx)P 
which, specified for the above parameters, yields 

„ fd — 3^^__1 ddd 



Fs{ai, 02, "2, A, /52, /?3; 7i, 7i, 7i; a;, ^) 



(92) 



2 '^'^'^'^'2'2'2'2''^'^'^. 





i) 




r 






(1) 







(1 — X + tx) 



2' 2' 



By means of 



Fi (a, 6, 6', 6 + b'; w, z) = (1 - ^F, 
Fi in the integrand reduces to 

1 



a,b; w — z 
b + b'; 



1 3 

^1 ( 1,1, 2' 2'*^^'^^ 



l-tz 



2F1 



1,1; tiy-z) 
3 . 

2 ' 



l-tz 



arcsm 



t(y-z) 

l-tz 



^t{y~z){l-ty) 



(93) 



(94) 



(95) 



and finally we obtain 



(i — 3 1 d d d ^ 

,1,1; 1,1, 2: 



4-5 
2 



r(f) (y-^)-^ .i arcsinv/B(l-t) 



dt. 



(96) 



Secondly the contribution to from 2F1 in (ffH|l can be evaluated in the same manner as 



in Sect. 5, where I^f^ was calculated from J2"''. The result reads 



Ad) 



/ diXijki\ ( dkX 



A 



i'jkl 
ijkl 



r r (t2) 



2 

7_ 



i?^ ( 1 1 ^~ 1 . IjL liL 

\ 2 2 2 l^ijkl ^ijk ^ 



(97) 



In complete analogy to (j74|) we now write the result for the 4-point function as 

&4 



r(2-f) 



+ (f^ijki diXijki + 4>iijk dkXijki + djXijki + 4>jkii diXijki, (98) 
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where 



v/2 



Hjkl 



1 ( dkXi 



\/~9ijk 



Kjk J 



1 ij— 3 . 

-L, 2 1 ' 

d-2 . 

2 ' f^ijkl 

d-2 ( 



1-7^ ^n 



\ 2 2 2 T^ijkl fijk, 

rf-2 ^ ^ ^ 



d-2 




'2 '2 
^ijkl ^ijk ^ij , 



k,i,j} + {ij,k 



j,k,i} 



(99) 



The first term on the r.h.s. of comes from 63 of (fTij) . but is supposed to be if the 
3— point function is evaluated in a kinematical domain, where no boundary term occurs. 
The boundary term 64 we give for completeness: 



3 

in 2 



r 


it 




{d-2)T{ 





d-3 
2 



9ijkl ^ijkl 



(100) 



with the same reasoning as for (j75j) . So far this result demonstrates our general method 
also for the 4-point function. 

To give an explicit and complete example we consider now the case with the following 
kinematical variables 



p\ = m^, pI = s, pI = m^, P1P2 = P2P3 = ^s, pips = 



mi 



mg = 0, m: 



(101) 



which corresponds to the (scalar) box diagram with two photons in the s-channel, occurring 
in Bhabha scattering. Apart from (jHTj) we also need the 3-point function Jp^^ with m\ = 
m\ = m^, m\ = and Pi = P2 = ^'^1 which is 



J 



id) 



r(2-i) 



2 m 



6-d 



1,3-f; q 
3 . 



2 ' 



4m^ 



(102) 



With these two 3-point functions we can set up the difference equation for the 4-point 
function with kinematics given in (jlOU) . the solution of which finally reads (the indices of 
/nil being the powers of the scalar propagators, see also (HSI)) 



Ad) 

'1111 
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t{s - 4m2) I 2 1 



'd-3 3 d-2 

-L 1 -L 1 ~r 



2 ' ' '2' 2 's-4m2^ 

d-3. d-3 



+ 



Am' 



d-4 



(rf - 3)t(s - 4m2) 



2 ■ 2 
d-1 . d-2 
2 ■ 2 



1; 1; 



m 1 



with 2; 



4tt 



tUni^-s) 



to 

geometric functions are 



and s, t, u the usual Mandelstam variables. Here the generalized hyper- 



'd-3 1 d-1 
^~'^'2'^~' 



00 00 (^) (I) 



V ^ / r+s 



(104) 



^\ 2 ' ' ' 2 2 



and the Kampe de Feriet function 



EE 

r=0 s=0 



r+s 



d-2 
2 



(105) 



pl;2;l 
-^i;i;0 



2 ■ 2 ' ' ' 
d-1 d-2 



x,y 



EE 

r=0 s=0 



d-3 
2 

d-1 
2 



d-3 
2 



r+s 



d-2 
2 



{x,y). (106) 



64 is the boundary term, which we have to determine in the following. First of all we 
calculate the derivatives dk^i- 



di/^A = <93A4 = 2st{Am^ - s), 



and 



92A4 = 54A4 



-2st^. 



(107) 



(108) 



These two sets of derivatives obviously have opposite sign for physical values of the kine- 
matical variables s and t (and also under their exchange for the crossed diagram). Thus, 
according to the discussion after (jSH), we conclude that the boundary term is zero. 

Nevertheless it is of interest to investigate an alternative approach to determine the 
boundary term: instead of the asymptotic method described in Sect. 3 we set up a differ- 
ential equation for 64 with respect to . An expression for the derivative of with 
respect to can be obtained by differentiating the parametric representation of / 



id) 
1111 



Ad) 
-'1111 



1 

'jd/2 



ndcy. -j 
—I exp 



■ Q . \ ^ 2 

« n - « E «fc"^fc 
. ^ k=i 



(109) 



where 



2 2 2 2 2 2 

Q = aia2Pi2 + aiftaPis + «i«4Pi + a2a3P23 + «2a4P2 + «3«4P3, 

D = tti + 0:2 + as + a^. (110) 
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Substituting the kinematical variables (jlUl|) gives 



Q — D ^ dkiTik = Oio^t + a2a4S — (02 + a^] 



2 2 
m . 



(Ill) 



k=l 



Differentiating p09|) we obtain 



arid) 

^-'1111 _ r, T{d+2) r) T(d+2) r) j{d+2) 
~ '^-'1311 ^-'1212 '^-'1113 • 



(112) 



Using recursions with respect to the indices and dimension as described in j7], the integrals 
on the right hand side can be reduced to master integrals with the result 



1111 



dm? u{Am?' — s 



\ ^{2[(4-rf)t-«]/Si-8(rf-3)/ 



(d) 



+ ?^lf{,n^) - - (d - 4) (^i±|^ Jg«'^ , (113) 



where Ip^ and Iq' are 2-point functions 



rid) 



Ad) _ Ad) 

-'g — -'2 



^ r(2-f)r(f-i) 



mi=0,m2=0 ('_j52-)(2-|) 2d-3 p (^^—^^ 



(d) _ Ad) 
F 



TW _ T 
J F — -'2 



[m 



mi=m,m2=m 



1,2-1; p2 



^' " 2 ' 
3 . 
2 ' 



4m^ 



(114) 



(115) 



Differentiating p03p w.r.t. and equating the two expressions for the derivative , 
substituting (|1U3|) into the r.h.s. of (jll3p yields a differential equation for 64 . To derive this 
equation we used 



F2 



2{x — y)x 



2{x — y)x 



F2 

dFi 
dx 

dFi 
dy 



'd-2 1 d ^ 
-^,l,-,-;x,y^ 

'd-3 3 d-2 
' ' ' o' 



^ 2 ' ' '2' 2 
= [x{2-d) + y{d-3)]Fi-y{d-3)Fiy 

+ x{l~y){d-2)F,[-,l,-,-;x,y] , 



'd-2 Id ^ 
-^,0,-,-;x,y 



yF, + y{d-3)F, 



r ' 2' 2' 

'rf-2 1 d 
,0, 



2 '"'2'2''^'^. 



Fi(a,0,/5',7; x, y) 
Fi{a,/3,/3\a;x,y) 



y{l - x){d - 2)Fi ^^,1,^, 

2Fi(a,/3',7,?/), 
1 



{l-xYil-yf 
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(116) 



For the Kampe de Feriet series the following relation was used 

ocbix v) 

2(1/ - x)x — ^ — = [y - (d - 3){y - x)](j){x, y) 

-x{d - 3)2Fi h , x\+{d- A)y 2F1 h , y j . (117) 

For the F2 function we used the relation: 

—8m stu— — = —8m (4m [t + u) — us)u— — 
ox dy 

+ t(4m^ - s)[(4m^(t + u) - us){d - 4) - 4m^M]F2 
+ tu{s - 8m2)(s - 4m^) 2^1 (^1, -m^z^ 

-(rf-4)t^(.-4mY2F,(l,3-^,|£,) (118) 
Finally the following differential equation was obtained for 64: 

dm? u{s — Am?) ^ ^ ^ 

As initial value we chose hi{rn? = 0) = 0, as we found explicitely and thus 64 = for also 
for w? 7^ 0. 

7 Conclusion 

For many two-loop problems in the electroweak theory it is necessary to calculate one-loop 
diagrams in arbitrary dimension, e.g. to obtain the e expansion when inserting counter- 
terms. For this problem we offer a general solution for 2—, 3— and 4— point functions for 
arbitrary masses and kinematics. The solution is obtained in terms of generalized hyper- 
geometric functions plus a 'boundary term'. To determine the latter, two different methods 
have been applied: an asymptotic expansion in the dimension d and alternatively a first 
order differential equation. Our results (for arbitrary kinematics) seem fairly lengthy, but 
they are finally expressed in terms of functions 2Fi,Fi and Fs, which are quite accessible 
analytically as well as numerically. In particular we point out the simple integral represen- 
tations ()78|1 and ()9fij) . which can easily be evaluated for arbitrary kinematics, and thus will 
become of particular importance for the evaluation of 5— point functions - as mentioned 
in the introduction. For diagrams occurring in Bhabha scattering we have given explicit 
results for arbitrary dimension of the needed one-loop diagrams. The obtained results will 
serve in further calculations as useful tools. 
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